Dispelling the curse of the neutron skin in atomic parity violation 
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We perform calculations for the neutron skin of nuclei and its contribution to atomic parity non- 
conservation (PNC) in many isotopes of Cs, Ba, Sm, Dy, Yb, Tl, Pb, Bi, Fr, Ra. Three problems are 
addressed: i) Neutron-skin induced errors to single-isotope PNC, ii) Possibility to measure neutron 
skin using atomic PNC, iii) Neutron-skin induced errors for ratios of PNC effects in different isotopes. 
In the latter case the correlations in the neutron skin values for different isotopes lead to cancelations 
of the errors; this makes the isotopic ratio method a competitive tool in a search for new physics 
beyond the standard model. 
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Atomic parity non-conservation (PNC) provides pow- 
erful constraints on extensions to the standard model of 
elementary particles in the low-energy electroweak sector. 
In such measurements one determines a parity- violating 
signal -EpNCi related to the quantity of interest, the weak 
charge, Qw, as i?pNC = ^pnc Qw- The coefficient fcpNc 
comes from atomic calculations. Considering challenges 
faced by such calculations, an alternative approach was 
proposed by Dzuba et al. [l| . The idea was to form a ratio 
TZ of the PNC amplitudes for two isotopes of the same ele- 
ment. Since the factor fcpNc remains the same, it cancels 
out in the ratio. However, Fortson et al. pointed out 
a conceptual limitation to this approach - an enhanced 
sensitivity of possible constraints on "new physics" to un- 
certainties in the neutron distributions. This problem is 
usually referred to as the problem of the neutron "skin" . 
Almost for two decades this problem has persisted. Here 
we show that the neutron skins in different isotopes are 
correlated; this leads to a substantial cancelation in the 
neutron skin induced uncertainties in the PNC ratios. 
The use of modern experimental data and nuclear cal- 
culations makes the isotopic ratio method a competitive 
tool in search for new physics beyond the standard model. 
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The neutron skin ARnp is defined as a difference be- 
tween the root-mean-square radii Rn and Rp of neutron 
and proton distributions. Even in interpreting the most 
accurate to date single-isotope measurement in Cs 3], 
this was a point of concern, as the induced uncertainty 
was comparable to experimental error bar for the PNC 
amplitude d, Q. The question was addressed in Ref. 
where empirical antiprotonic-atom data fit for the neu- 
tron skin was used [Vj, and the associated uncertainty 
in the "skin" contribution to Ep mc was substantially re- 
duced. At the same time this question is yet to be settled 
for on-going PNC experiments with unstable analogs of 
Cs: Fr 8] and Ra+ 9]. In this letter we use results of 
recent advances in our theoretical and experimental un- 
derstanding of neutron skins to address the important 
questions that can be answered from atomic parity vio- 
lation experiments. 
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FIG. 1: The neutron skin AR„p for nuclei above ^°°Sn are 
shown in the bottom panel (all Z, even A'^ values). The nuclei 
for a given isotope are connected by lines. The filled circles 
are those for the nuclei of interest for atomic parity violation. 



Nuclear- structure calculations- The nuclear structure 
models used in the present calculations were constrained 
by a number of measured observables, primarily by the 
anti-protonic data. Until recently the anti-protonic atom 
data were analyzed in terms of a simple Fermi shape for 
the matter density. The anti-protonic data are sensitive 
to the very low-density part of the matter distribution, 
for •^"^Pb, for the density near i? = 10 fm [13] • One can- 
not justify using the simple Fermi shape for the density 
distribution at such large radii. 

Recently the anti-protonic atom data have been ana- 
lyzed in terms of nuclear EDF models 10] including those 
for the Skyrme and Relativistic Hartree formulations. At 
large radii the matter density is dominated by the neu- 
tron density (due to the Coulomb barrier for protons). 
But the connection between the asymptotic density and 
the rms radius within the energy-density functional mod- 
els depends on other features of models such as the nu- 
clear matter incompressibility K . The charge density ob- 
tained from electron scattering is best reproduced with 
EDF models with if w 200 - 230 MeV. This eliminates 
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many of the parameter sets which give higher values of K. 
Within the Skyrme models K is closely controlled by the 
power of the density-dependent potential, p", with a=l 
with K « 330 down to a=l/6 with K Ri 200. The anti- 
protonic atom data were compared to results of three 
new Skyrme forces called Skxsxx with a=l/6 and 
values of xx=15, 20 and 25 representing the neutron skin 
of ■^°^Pb in units of 10~^fm. The result of the analy- 
sis was Ai?„p = 0.20(±0.04)(±0.05) fm, where ±0.04 fm 
is experimental error from the anti-protonic line width, 
and where ±0.05 fm is the theoretical error suggested 
from the comparison of the theoretical and experimen- 
tal charge densities at large radii. The neutron skin can 
also be constrained by the properties of the pygmy dipole 
resonance in neutron- rich nuclei Data for ^^^Sn sug- 
gest a value of Ai?„p = 0.24(4) /m This is consistent 
with the value of 0.27(5) fm obtained for ^'^^Sn from the 
Skxs20(5) interactions. 

The neutron skins for nuclei above ^'^''Sn obtained with 
Skxs20 are shown in Fig. 1. These are obtained in a 
spherical basis with nuclcons allowed to occupy the lowest 
energy orbitals - the calculation does not have deforma- 
tion or pairing. The irregularities between the neutron 
magic numbers 82 and 126 are due to the filling of the 
proton hii/2 orbital and the neutron ii3/2 orbital. How- 
ever, since the orbitals in the major shells (g7/2, d5/2, 
d3/2, Si/2, hn/2) for protons and (h9/2, f7/2, P3/2, 
Pi/2; 113/2) for neutrons are closely spaced, deformation 
and pairing will average out these irregularities. Thus, 
in the discussion of the values and their errors we will 
use average values obtained from a smoothed variation 
in the neutron-skin with neutron number. The values 
can be improved with calculations that include pairing 
and deformation, but our conclusions should not change. 
The error bars for the neutron skin were obtained from 
[AR„p(Skxs25) - AR„p(Skxsl5)]/2. 

Atomic PNC and neutron skin ~ The PNC observ- 
ables depend on matrix elements of weak interaction. As 
demonstrated in Ref. [13], the matrix elements of i7w 
may be parameterized as 



The accuracy of the above formula is sufficient for the 
goals of the present work. 

Single-isotope measurements - The relative correction 
to the PNC amplitude due to the neutron skin reads 0] , 
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The computed corrections for all the isotopes are listed 
in TableU In particular, for ^^"^Cs the relative correction 
is -0.0023(5), it is consistent with the value of -0.0019(8) 
from Ref. 6], which was based on the semi-empirical fit 
of antiprotonic-atom data (Ai?„p = 0.13(4) fm). As 
we progress to heavier elements, the correction grows as 
Z^, reaching 0.6% for Fr and Ra"*". As an example, for 
^^■^Fr, the correction reads —0.0063(16). The error bar 
implies that at the present level of knowing the neutron 
skin, it contributes to the uncertainty in the extraction 
of new physics from the Fr experiment at the 0.1-0.2% 
level. 

Detecting neutron skin in atomic PNC — The ques- 
tion of determining neutron skin is of interest in its own 
right, for example, for equation of state for neutron stars. 
It is worth mentioning the proposed PREX experiment 
at JLAB [l3| on using a PNC asymmetry in elastic scat- 
tering of electrons from ^"^Pb to measure i?„ to a 1% (± 
0.05 fm) accuracy. There have been renewed attempts to 
obtain i?„ from hadronic scattering data [3, Ell ■ 

Considering this interest, we would like to see if the 
neutron skin can be extracted from atomic PNC mea- 
surements. From the preceding discussion, it is clear that 
for the single-isotope PNC the uncertainty of experiments 
and atomic calculations should be smaller than 0.2% (Cs, 
Ba"*") and 0.6% (Fr, Ra+). This seems to be a realistic 
goal [13,113. 

This problem may be also addressed in the isotopic 
chain experiments. Suppose the PNC amplitudes EpNc 
and E'pj^fj are measured for two isotopes of the same 
atom or ion with neutron numbers N and N' — + AA^, 
and the ratio is formed 
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where factor Cji depends on atomic wavefunctions, 7 = 
\/l — {c(Z)^ and Qw includes the dependence on nuclear 
distributions. 
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Here all quantities with primes are for the isotope with 
A^' neutrons. Focusing on the contribution of the neutron 
skin 



)w = ~Nqn + Zqp{l- Asm^0w) + AQ„ 
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The term AQ^ew characterizes "new physics" and 9w is 
the Weinberg angle. The quantities q„ and qp depend 
on the neutron and proton distributions convoluted with 

atomic wavefunctions: = 1 + /n (^^^ - In the "sharp 

edge" model of nuclear density distribution. 
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If we neglect the difference between i?„ and Rp, i.e., the 

neutron skin, then 7^ ^ 7^o = N/N' {Rp / R'pf^ ^'^ . Any 
deviation of TZ from TZq is a signature of the neutron skin. 
The figure of merit is 
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where /„ is given by Eq.(l3]). Or in terms of the neutron 
skin ATZn... « f [aZf ^ [Ai?;^ - Ai?„p] . 

There are two observations that can be made: (i) The 
isotopic-ratios are sensitive to the differential change in 
the skin thickness, i.e., the neutron-skin effects tend to 
cancel and (ii) the largest effect is attained for a pair of 
isotopes where the skin thicknesses differ the most. This 
condition is reached for a pair comprised of the lightest 
(neutron-depleted) and the heaviest (neutron-rich) iso- 
tope of the chain. 

For Cs, Ba, and Dy A7^„., w 0.001, A7^„.s(Yb, Sm) w 
0.002, and A7^„.^(Fr, Ra) « 0.003. Fr and Ra are the 
extreme cases, as the skin thicknesses for the lightest and 
the heaviest isotopes differ by a factor of two. 

Isotopic ratios: neutron skin vs. "new physics" - Sup- 
pose we form the ratio of measured PNC amplitudes for 
two isotopes of the same atom or ion. Can we constrain 
new physics beyond the Standard Model from analyz- 
ing the ratios? The previous studies have answered this 
question negatively, as the uncertainties in the neutron 
skin masked the new physics contributions. Below, in 
light of our nuclear calculations, we re- visit this question. 
The discussion follows previous analysis There are 
two new points: (i) we use more finely-tuned calcula- 
tions of the skin (TableJT] in lieu of the empirical fit of 
antiprotonic-atom data [7|) and (ii) we take into account 
that the errors in Ai?„p for two isotopes are correlated. 
Both these factors allow us to argue that, by contrast 
to the previous studies, the isotopic ratios can provide 
competitive constraints on the new physics. 
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FIG. 2: Neutron skin vs. "new physics". The neutron-skin- 
induced uncertainties for isotopic chains are compared with 
the constraints from parity-violating electron scattering. 

The term AQnew in Eq. ([2|) characterizes "new 
physics" at the tree level. Following we represent it 
as a combination of new-physics couplings to protons and 
neutrons, AQncw = Z hp + N hn- Various elementary- 
particle scenarios for these interactions were reviewed in 
Ref. [2(|. Then Qw = N ho + Z hp + N hn, where ho 
comes from the Standard Model (this also includes the 
nuclear corrections). Unlike in the single-isotope mea- 
surements (sensitive mainly to hn) in the isotopic ratio 



method, the sensitivity to new physics comes predomi- 
nantly due to hp. The sensitivity may be parameterized 
as M\ 
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In the absence of new couplings J- = 0. The smaller JF, 
the tighter the constraints on hp are. For a given chain, 
it is beneficial to work with the largest possible neutron 
spread AiV, i.e., forming the pairs from the lightest and 
the heaviest elements of the chain. 

The constraints on hp, Eq. ([7]), are affected by (i) the 
experimental errors, STZ^xp and (ii) uncertainties in TZo 
which are induced by insufhcient knowledge of nuclear 
distributions. Explicitly, 
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At this point we focus on the error in the difference 
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An important point is that the neutron-skin errors in 
Eq.® for two isotopes are correlated. From numerical 
experimentation, we find that for a given isotopic chain a 
variation of nuclear interaction parameters induces simi- 
lar change in the ratio Rp for all the isotopes. Indeed, 
in Table U the boundaries of error bars in the values of 
neutron skin correspond to the same values of nuclear 
parameters. We see that the errors tend to cancel each 
other. This is to be contrasted with the previous anal- 
ysis of Ref. [1^: there the errors from individual iso- 
topes were treated independently, i.e., errors were added 
in quadratures. We find that our new "correlated" treat- 
ment reduces the error bars by four to ten times. 

We calculated 5T from a smoothed variation in the 
neutron-skin with neutron number over pairs of the ad- 
jacent isotopes (AA^=2). The results are compiled in 
Table [II] and also in Fig. [2] These errors are compared 
to the new-physics couplings to protons, hp, which are 
directly probed by the parity-violating electron scatter- 
ing (PVES). The PVES experiments were recently ana- 
lyzed in [21]. The resulting weak charge of the proton is 
I22II. The error bar determines the 



= 0.058 ± 0.023 _ 
upper bound on T shown in Fig. (5] It is clear that all 
isotopic-chain determinations are competitive to bounds 
derived from PVES. For example, measurements with 
isotopes of Cs, Ba and Dy would be an order of magni- 
tude more sensitive to the new physics. 
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TABLE I: Computed neutron skins Ai?„p and proton distri- 
bution rms radii. The relative contributions of the neutron 
skin to atomic PNC amplitudes are listed in the last column. 



A 


ARrip, fm 


Rp, fm 




Cs {Z 


= 55) 






129 


0.120(32) 


4.705 


-0.0018(5) 


131 


0.139(34) 


4.714 


-0.0020(5) 


133 


0.158(37) 


4.723 


-0.0023(5) 


135 


0.176(40) 


4.732 


-0.0026(6) 


137 


0.193(42) 


4.742 


-0.0028(6) 


Ba {Z 


= 56) 






130 


0.104(29) 


4.723 


-0.0016(4) 


132 


0.124(32) 


4.731 


-0.0019(5) 


134 


0.143(34) 


4.740 


-0.0022(5) 


136 


0.161(37) 


4.749 


-0.0024(6) 


138 


0.179(40) 


4.759 


-0.0027(6) 


Sm {Z 


= 62) 






144 


0.098(27) 


4.894 


-0.0018(5) 


146 


0.122(32) 


4.893 


-0.0022(6) 


148 


0.144(35) 


4.903 


-0.0026(6) 


150 


0.166(39) 


4.913 


-0.0030(7) 


152 


0.187(43) 


4.924 


-0.0033(8) 


154 


0.219(48) 


4.934 


-0.0039(9) 


Dy (Z 


= 66) 






156 


0.135(35) 


4.997 


-0.0027(7) 


158 


0.150(37) 


5.018 


-0.0030(7) 


160 


0.164(38) 


5.039 


-0.0032(8) 


162 


0.178(40) 


5.061 


-0.0035(8) 


164 


0.191(42) 


5.082 


-0.0037(8) 


Yb {Z 


= 70) 






168 


0.141(35) 


5.143 


-0.0031(8) 


170 


0.153(38) 


5.163 


-0.0033(8) 


172 


0.174(40) 


5.171 


-0.0038(9) 


174 


0.202(51) 


5.173 


-0.0044(11) 


176 


0.215(67) 


5.193 


-0.0046(14) 


Tl {Z -- 


= 81) 






203 


0.179(45) 


5.422 


-0.0049(12) 


205 


0.192(48) 


5.434 


-0.0053(13) 


Pb {Z 


= 82) 






204 


0.172(44) 


5.430 


-0.0049(12) 


206 


0.184(46) 


5.442 


-0.0052(13) 


208 


0.200(50) 


5.450 


-0.0056(14) 


Bi (Z = 


= 83) 






209 


0.189(48) 


5.468 


-0.0054(14) 


Fr {Z = 


= 87) 






209 


0.121(36) 


5.518 


-0.0038(11) 


211 


0.132(38) 


5.529 


-0.0041(12) 


213 


0.146(42) 


5.536 


-0.0046(13) 


215 


0.161(44) 


5.546 


-0.0050(14) 


217 


0.176(47) 


5.555 


-0.0055(15) 


219 


0.191(50) 


5.565 


-0.0059(16) 


221 


0.206(53) 


5.574 


-0.0064(16) 


Ra {Z 


= 88) 






210 


0.111(34) 


5.535 


-0.0035(11) 


212 


0.123(37) 


5.546 


-0.0039(12) 


214 


0.136(40) 


5.553 


-0.0043(13) 


216 


0.151(43) 


5.563 


-0.0048(14) 


218 


0.166(46) 


5.572 


-0.0053(15) 


220 


0.181(49) 


5.581 


-0.0057(16) 


222 


0.195(52) 


5.591 


-0.0062(16) 
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TABLE II: Contribution of nuclear-structure uncertainty to 
a constraint on "new physics" obtained from a smoothed 
variation in the neutron-skin with neutron number. 



Atom 


Alas 


s imuibeivs .1 


()J^ X 10- 


Cs (Z = 55) 


129 


137 


2.1 


Ba {Z = 56) 


130 


138 


2.3 


Sm {Z = 62) 


144 


154 


4.2 


Dy {Z = 66) 


156 


164 


2.7 


Yb (Z = 70) 


168 


176 


10.2 


Tl (Z = 81) 


203 


205 


7.2 


Pb {Z = 82) 


204 


208 


7.7 


Fr {Z = 87) 


209 


221 


8.8 


Ra {Z = 88) 


210 


222 


8.9 



